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We discuss a new scenario for early cosmology when the inflationary de Sitter phase emerges 
dynamically. This genuine quantum effect occurs as a result of dynamics of the topologically non¬ 
trivial sectors in a strongly coupled QCD- like gauge theory in an expanding universe. We test 
these ideas by explicit computations in hyperbolic space x §Ui. We argue that the key element 
for this idea to work is the presence of nontrivial holonomy computed along S^-i. The effect is 
non-local in nature, non-analytical in coupling constant and can not be described in terms of any 
local propagating degree of freedom such as scalar inflaton field $(*). We discuss some profound 
phenomenological consequences of this scenario for inflationary cosmology. We also suggest to test 
these Ideas in a tabletop experiment by measuring some specific corrections to the Casimir pressure 
in the Maxwell theory formulated on a topologically nontrivial manifold. 


I. INTRODUCTION. MOTIVATION. 


The main motivation for the present studies is the pro¬ 
posal that inflationary de Sitter phase [1-3] may be dy¬ 
namically generated as a result of presence of the topolog¬ 
ically nontrivial sectors in expanding universe. Inflaton 
in this framework [4, 5] is an auxiliary topological non¬ 
propagating field with no canonical kinetic term, similar 
to known topologically ordered phases in condensed mat¬ 
ter systems. This auxiliary field effectively describes the 
dynamics of the topological sectors | fc) in a gauge theory 
(coined as QCD in [4, 5]) in expanding Universe. 

This picture should be contrasted with conventional 
proposals reviewed in [2, 3] when the de Sitter behaviour 
is achieved in quantum field theory (QFT) by assuming 
the existence of a new scalar local field <I’(x) with a non¬ 
vanishing potential energy density U(4>). The shape of 
this potential energy can be adjusted in a such a way 
that the contribution to energy density e and pressure p 
is in agreement with observations. In different words, the 
scale parameter a(t) and the equation of state during the 
inflation take the following approximate form. 


a(t) ~ exp{Ht), e « —p. (1) 

The key ingredient of the proposal [4, 5] is a conjecture 
that the vacuum energy in context of the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) Universe has the 
following expansion at small H <C Aq^q 
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when the first non-vanishing term is linear ~ H, rather 
than (commonly accepted) quadratic ^ in the Hubble 
constant. If this conjecture turns out to be correct, than 
the Friedman equation assumes the form 
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which automatically leads to a non-trivial solution with 
constant Hq, and as a consequence, to a desired de Sitter 
behaviour (1). 


There are two critical elements in writing equation 
(3). First one, as we already mentioned, is related to 
the expansion (2), see few comments on this conjecture 
below. The second key element is a paradigm that the 
relevant definition of the energy in an expanding back¬ 
ground which enters the Friedman equation is the differ¬ 
ence AE{H) = [E{H) — AMink], similar to the computa¬ 
tion of the Casimir pressure when the observable energy 
is the difference similar to AE. This element in our anal¬ 
ysis is not a new proposal identifying AE with gravitat¬ 
ing energy from the Friedman equation. In fact, in the 
present context such a definition for the vacuum energy 
was advocated long ago in 1967 by Zeldovich [6] for the 
first time. Later on such definition for the relevant energy 
AE = (ApLRW ~ ^'Mink) which enters the Friedman equa¬ 
tions has been advocated from different perspectives in a 
number of papers, see e.g. relatively recent works [7-11], 
see also review article [12] with large number of references 
on original papers. Essentially, this prescription implies 
that AE may only depend on properties of the external 
gravitational background, while the conventional contri¬ 
butions computed in Minkowski flat space-time (such as 
the QCD vacuum energy or the Higgs potential in elec- 
troweak theory) are automatically subtracted by this pre¬ 
scription^ . We shall not elaborate on a number of subtle 
points related to this prescription in the present work 
referring to the original papers and review article [12]. 

The main topic of the preset paper is analysis of an- 


^ A somewhat similar, but not identically the same subtraction 
procedure has been suggested recently in refs. [13, 14], the so- 
called “vacuum energy sequestering” proposal. The prescription 
[13, 14] is also inherently non-local, similar to the crucial role of 
non-locality in our framework realized in terms of the holonomy 
(5). In fact, our computation of the vacuum energy (2) as dis¬ 
cussed below, is based on evaluation of the holonomy (5) along 
the entire history of the universe, which resembles in spirit the 
computations of the so-called “historic averages” in refs. [13, 14]. 
Furthermore, we have to keep the volume of the system to be 
finite in the computations for the infrared regularization of the 
theory. It is akin to that from refs. [13, 14] where the finite volume 
is also required property for the consistency of the procedure. 
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other critical element, briefly mentioned above, and lead¬ 
ing to (3). There is well known, conventional and gener¬ 
ally accepted argument which suggests that the expan¬ 
sion (2) starts with a quadratic ~ H^, rather than the 
liner ~ H term. The argument is based on fundamental 
principles of locality and general covariance, see original 
papers [15, 16], recent review [12], and some comments 
[4] with pros and cons of these arguments. Indeed, the 
curvature scalar R for FLRW Universe is quadratic in H, 

^ + 6H, (4) 

when H ~ see [12]. Therefore, if the in¬ 

frared (IR) behaviour of the system is entirely determined 
by the local characteristics, such as curvature scalar R 
and/or higher order derivative terms i?^, than 

the corrections to the energy (2) indeed must be pro¬ 
portional to even powers as correctly argued in 

[12, 15, 16]. 

However, the main essence of the proposal [4, 5] is pre¬ 
cisely the observation that the conventional assumption 
on locality might be badly violated in strongly coupled 
gauge theories. The basic reason for such violation is well 
known and well-understood, at least in Minkowski space- 
time. The energy (2) is generated due to the tunnelling 
events between |fc) topological sectors, which formulated 
in terms of inherently non-local large gauge transfor¬ 
mation operator T. Furthermore, this energy has non- 
dispersive nature, i.e. it can not be formulated in terms 
of any local propagating degrees of freedom^. Transi¬ 
tion from Minkowski space-time to time dependent back¬ 
ground (2) obviously will not modify the nature and ori¬ 
gin of this type of energy. Rather, a transition to FLWR 
Universe introduces some background- dependent correc¬ 
tions to the same type of energy (2), which was coined 
as “strange energy” in [4, 5] due to its unconventional 
origin as mentioned above. 

One should comment here that this feature of non¬ 
locality when the system is not completely characterized 
by a local physics is very similar to the well known prop¬ 
erty in topologically ordered phases in condensed matter 
physics wherein an expectation value of a local operator 
does not fully characterize the ground state of the sys¬ 
tem. Instead, one should use some non-local variables for 
proper characterization of the system. 

The main subject of the present work is to elaborate 
and clarify a number of non-trivial questions related to 
the non-locality in QFT and generation of the linear ^ H 
term in (2) in some simplified models^. The basic point 
of our discussions is that a gauge QFT (when one should 


|R|-6(^^ + 


sum over all topological sectors |fc) in the definition of 
the partition function) is not fully described by the lo¬ 
cal characteristics, such as curvature (4). In particular, 
the linear dependence on the background may enter (2) 
through other characteristics of the system such as holon- 
omy 


C/(x) = R exp 



dxiAi{xi,yi) 


(5) 


which is gauge invariant but non-local object as it de¬ 
pends on the boundary conditions. We shall argue below 
that precisely the non-trivial holonomy in gauge theories 
plays a key role in the mechanism which could gener¬ 
ate the “strange energy” (2). It is very hard technical 
problem to compute the non-perturbative energy (2) in 
a time- dependent background characterized by parame¬ 
ter H, see footnote 3 for clarification. However, one can 
simplify the problem by considering the sensitivity of a 
gauge system to some external dimensional parameters 
characterizing the gravitational background, such as n, 
see definition below. This parameter plays a role sim¬ 
ilar to the Hubble constant H in FLRW Universe (2). 
Our goal is to study the dependence of the “strange en¬ 
ergy” (2) as a function of n in the limit of small k —?► 0 
in some simple settings where such computations can be 
performed. 

The basic idea is as follows. We would like to consider 
hyperbolic space with the constant negative curvature 
As we discuss below, there is a conformal equiv¬ 
alence between and x S/_i where 

denotes the circle of radius . The holonomy (5) is 
computed precisely along a closed loop Our goal 

is to study the first nontrivial correction ~ k to the non- 
perturbative energy (2) in the limit of small k —>■ 0 cor¬ 
responding to smooth transition to Euclidean space R^. 
We should recover the Euclidean results when k identi¬ 
cally vanishes. 

The key observation is that the topological configu¬ 
rations with non-trivial holonomy (5) produce a finite 
contribution to the energy density (2) with corrections 
being linearly proportional to k. Such effects can not be 
expressed in terms of any local operators such as cur¬ 
vature (4). Rather, it is generated due to topological 
vacuum configurations with nontrivial holonomy (5), not 
expressible in terms of local observables. This is precisely 
the reason why the generic arguments [12, 15, 16] based 
on locality simply do not apply here. 

Therefore, our computations of the linear correction 
~ K in the vacuum energy density using simplified model 


^ This energy can be expressed in terms of the contact term in 
the topological susceptibility, determined by the IR physics and 
boundary conditions. The corresponding physics has been well 
understood using the lattice numerical simulations in strong cou¬ 
pling regime, see [4] for references and details. 

® Here and in what follows we use term “linear in iJ correction” 


as a generic feature to distinguish a nontrivial background from 
the trivial Euclidean space. This dimensional parameter should 
not be literally identified with the Hubble constant. Rather, it 
could be any other dimensional parameter which characterizes 
the system, such as the size of torus related to the nontrivial 
holonomy (5) with p ~ H~^. 
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with X S^_i background can be thought as a strong 
argument supporting our conjecture on the liner correc¬ 
tion ~ H in generic FLRW Universe (2)- in both cases 
the linear correction is not associated with the local cur¬ 
vature operator (4). 

One may wonder how a bulk property (such as vacuum 
energy density) in a gapped theory could be ever sensi¬ 
tive to such global characteristics as radius of the circle of 
The answer lies not in the local, but global proper¬ 
ties of the space. Imagine that we study the Aharonov- 
Casher effect. We insert an external charge into a su¬ 
perconductor in which the electric field is exponentially 
suppressed ~ exp(—r/A) with A being the penetration 
depth. Nevertheless, a neutral magnetic fluxon will be 
still sensitive to an inserted external charge at arbitrary 
large distances in spite of the screening of the physical 
field (which is equivalent to the presence of a gap in our 
system). This genuine quantum effect is purely topologi¬ 
cal and non-local in nature and can be explained in terms 
of the dynamics of the gauge sectors which are respon¬ 
sible for the long range dynamics. Imagine now that we 
study the same effect but in a different background. The 
corresponding topological sectors will be modified due to 
the variation of the external background. However, this 
modification can not be described in terms of any local 
dynamical fields, as there are no any propagating long 
range helds in the system since the physical electric field 
is screened. For this simplified example, the dynamics of 
the “strange energy” as a function of k is determined by 
the modifications of topological sectors when the back¬ 
ground varies. The effect is obviously non-local in nature 
as the Aharonov-Casher effect itself is a non-local phe¬ 
nomenon. 

The readers interested in the cosmological conse¬ 
quences, rather than in technical computational details 
may directly jump to section III D where we list the main 
results of this framework. Section IV is our conclusion 
where we suggest to test some of the ideas presented in 
this work in a tabletop experiment by measuring some 
specific corrections to the Casimir vacuum energy in the 
Maxwell theory formulated on 4-torus with the nontrivial 
holonomy. 


II. THE NATURE OF “STRANGE” ENERGY 
(2) IN EUCLIDEAN SPACE 

The main goal of this section is to review a number of 
crucial elements relevant for our studies. We start in sec¬ 
tion IIA with explanation of a highly nontrivial nature of 
the “strange energy” (2) in the Euclidean space time. We 
continue in section IIB, by clarifying the crucial role of 
the holonomy (5) in generating such type of energy. We 
present few technical results in sections IIC, IID. Finally, 
in section IIE we elaborate on non-local features of rel¬ 
evant vacuum conhgurations saturating the “strange en¬ 
ergy” in cosmological context. The corresponding anal¬ 
ysis will play an important role in our main section III 


when we study the hyperbolic space x S).-i and ana¬ 
lyze the properties of the “strange” energy as a function 
of K at small k —>■ 0. 


A. The topological susceptibility and contact term 


We start our short overview on the “strange” nature 
of the vacuum energy (2) by reviewing a naively unre¬ 
lated topic- the formulation and resolution of the so- 
called U{1)a problem in strongly coupled QCD [17-19]. 
We introduce the topological susceptibility % which is ul¬ 
timately related to the vacuum energy E^s^{9 = 0) as 
follows^ 

= lim / d^xe*'=“(T{g(a;),q(0)})(6) 
e=o '=^0 J 

where 9 parameter enters the Lagrangian along with 
topological density operator q{x) = and 

Eva,c{d) is the “strange” vacuum energy density which 
represents the first term in expansion (2) corresponding 
to the flat space-time background. This 9- dependent 
portion of the vacuum energy (computed at 9 = 0) has 
a number of unusual properties as we review below. The 
corresponding properties are easier to explain in terms of 
the correlation function (6), rather than in terms of the 
vacuum energy Fva.c{9 = 0) itself. The relation between 
the two is given by eq. (6). 

First of all, the topological susceptibility y does not 
vanish in spite of the fact that q = dfj,K^ is total diver¬ 
gence. This feature is very different from any conven¬ 
tional correlation functions which normally must vanish 
at zero momentum if the corresponding operator can be 
represented as total divergence. 

Secondly, any physical |n) state gives a negative con¬ 
tribution to this diagonal correlation function 


5^Evac(0) 

902 


Xdispersive ^ lim f {T{q{x), q{0)}) 

k^oj 

—kr — mi ^ mi 

n n 


(7) 


where is the mass of a physical jn) state, fc —> 0 is 
its momentum, and (Ojgjn) = c„ is its coupling to topo¬ 
logical density operator q{x). At the same time the res¬ 
olution of the U{1)a problem requires a positive sign for 
the topological susceptibility (6), see the original refer¬ 
ence [19] for a thorough discussion, 


Xnon-dispersive = lim [ xe"’^^{T{q{x), q{0)}) > 0. (8) 

k^Oj 

Therefore, there must be a contact contribution to X; 
which is not related to any propagating physical degrees 


^ We use the Euclidean notations where path integral computa¬ 
tions are normally performed. 
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of freedom, and it must have the “wrong” sign. The 
“wrong” sign in this paper implies a sign which is oppo¬ 
site to any contributions related to the physical propa¬ 
gating degrees of freedom (7). The “strange energy” in 
this paper implies the 9 dependent portion of the energy 
(2), (6) which can not be formulated in terms of con¬ 
ventional propagating degrees of freedom as it has pure 
non-dispersive nature according to eqs. (7), (8). 

In the framework [17] the contact term with “wrong” 
sign has been simply postulated, while in refs. [18, 19] 
the Veneziano ghost (with a “wrong” kinetic term) had 
been introduced into the theory to saturate the required 
property (8). 

Third, the contact term (8) has the structure x ~ 
J d‘^xS'^{x). The significance of this structure is that the 
gauge variant correlation function in momentum space 

lim j d^xe^^^{K^{x),K^{0)) - (9) 

develops a topologically protected “unphysical” pole 
which does not correspond to any propagating massless 
degrees of freedom, but nevertheless must be present in 
the system. Furthermore, the residue of this pole has the 
“wrong sign”. This “wrong sign” is due to the Veneziano 
ghost contribution saturating the non-dispersive term in 
gauge invariant correlation function (8), 

{q{x)q{0)) ~ {df,K^^{x),d^K''{0)) - S'^ix). (10) 

We conclude this review-type subsection with the follow¬ 
ing remark. The entire framework, including the singular 
behaviour of {q{x)q{0)) with the “wrong sign”, has been 
well confirmed by numerous lattice simulations in strong 
coupling regime, and it is accepted by the community as 
a standard resolution of the U{1)a problem. Further¬ 
more, it has been argued long ago in ref. [20] that the 
gauge theories may exhibit the “secret long range forces” 
expressed in terms of the correlation function (9) with 
topologically protected pole at fc = 0. 

Finally, in a weakly coupled gauge theory (the so-called 
“deformed QCD” model [21]) where all computations can 
be performed in theoretically controllable way one can 
explicitly test every single element of this entire frame¬ 
work, including the topologically protected pole (9), the 
contact term with “wrong sign”, etc, see ref. [22, 23] for 
the details. In particular, one can explicitly see that the 
Veneziano ghost is in fact an auxiliary topological field 
which saturates the vacuum energy and the topological 
susceptibility x- What is more important for the present 
studies is that one can explicitly see that the holonomy 
(5) plays a crucial role in generating the “strange” vac¬ 
uum energy and contact term in topological susceptibil¬ 
ity. 

While all these unusual features of the vacuum energy 
are well-known and well-supported by numerous lattice 
simulations (see e.g. [23] for a large number of references 
on original lattice results) the analytical understanding of 
these properties in strong coupling regime is still lacking. 
In next subsection we review some known results on this 


matter specifically emphasizing on role of the holonomy 
(5) in the analytical computations. Precisely a nontrivial 
holonomy (5) may play a crucial role in generating the 
linear correction ~ in eq. (2) as we argue in section 
HID. This is the key technical element which pinpoints 
the source of linear corrections ^ H not expressible in 
terms of any local operators such as curvature (4). 

In conclusion we should comment that the vacuum en¬ 
ergy in electroweak (EW) sector of the standard model 
is not sensitive to the the topological features of the EW 
gauge fields of the Standard Model (lF)f, A^) as these 
topological properties play no role in dynamics of these 
fields, in huge contrast with QCD. This is due to the 
fact that the mass gap of the non abelian EW gauge 
bosons is resulted from the Higgs mechanism, in contrast 
with QCD where the mass gap and the vacuum energy 
are dynamically generated by the topological fluctuations 
with nontrivial holonomy (5). Therefore, the subtraction 
procedure, formulated in Introduction, automatically re¬ 
moves all the vacuum energy related to the EW sector of 
the standard model without any corrections ^ H. The 
linear corrections are specific to the strongly coupled 
QCD with its topological features formulated in terms 
of the non-local operators, the holonomy (5) and large 
gauge transformations operator T as discussed above. 


B. The holonomy (5) and generation of the 
“strange” energy in Euclidean space. 


The key role in our discussions will play the behaviour 
of holonomy U (x) at spatial infinity, the Polyakov line. 


L = V exp 


dx 4 ^Ai{x 4 ^^ jxj —>■ oo) 


( 11 ) 


The operator TrL classifies the self-dual solutions which 
may contribute to the path integral at finite temperature 
T = including the low temperature limit T —^ 0. 
There is a well known generalization of the standard 
self-dual instantons to non-zero temperature, which cor¬ 
responds to the description on x geometry. This 
is so-called periodic instantons, or calorons[24] studied 
in details in [25]. These calorons have trivial holonomy, 
which implies that the TrL assumes values belonging to 
the group centre Z^r for the SU{N) gauge group. 

More general class of the self-dual solutions with non¬ 
trivial holonomy (11), the so-called KvBLL calorons were 
constructed more recently [26, 27]. In this case the holon¬ 
omy (11) in general, is not reduced to the group cen¬ 
tre TrL ^ Z^v. The fascinating feature of the KvBLL 
calorons is that they can be viewed as a set of N 
monopoles of N different types. Normally, one expects 
that monopoles come in iV — 1 different varieties carrying 
a unit magnetic charge from each of the U{1) factors of 
the gauge group left unbroken by vacuum ex¬ 

pectation value due to nontrivial holonomy (11). There 
is an additional, so called Kaluza- Klein (KK) monopole 
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which carries magnetic charges and instanton charge. All 
monopole’s charges are such that when complete set of 
different types of monopoles are present, the magnetic 
charges exactly cancel, and the configuration of N differ¬ 
ent monopoles carries a unit instanton charge. 

It has been known since [25] that the gauge configura¬ 
tions with non-trivial holonomy are strongly suppressed 
in the partition function. Therefore, naively KvBLL 
calorons can not produce a finite contribution to the par¬ 
tition function. However, this naive argument is based on 
consideration of the individual KvBLL caloron, or finite 
number of them. If one considers a grand canonical as¬ 
semble of these objects than their density is determined 
by the dynamics, and the old argument of ref. [25] breaks 
down. The corresponding objects in this case may in fact 
produce a finite contribution to the partition function. A 
self consistent computations in a weak coupling regime 
supporting this picture have been carried out in the so- 
called “deformed QCD” model [21]. One can explicitly 
see how N different types of monopoles with nontrivial 
holonomy (11) which carry fractional topological charge 
±l/N produce confinement, generate the “strange” vac¬ 
uum energy (6) and associated with this energy the topo¬ 
logical susceptibility (10) with known, but highly unusual 
properties reviewed above in section II, see [22, 23] for the 
technical details on these computations. 

In strong coupling regime we are interested in, the 
corresponding analytical computations have never been 
completed. There is a limited number of partial ana¬ 
lytical and numerical results [28-32] on computations of 
moduli space and one loop determinant, controlling the 
dynamics and interaction properties of the constituents 
in a large ensemble of KvBLL calorons. We review these 
basic technical results in next section IIC, as they will 
play an important role in our analysis below. 

While complete analytical solution in strong coupling 
regime is still lacking, nevertheless there is a number of 
hints supporting the basic picture that the KvBLL con¬ 
figurations with nontrivial holonomy (11) and represent¬ 
ing N different types of monopoles with fractional topo¬ 
logical charges ±l/iV saturate the “strange” vacuum en¬ 
ergy (6) and associated with this energy the topological 
susceptibility (10) in a very much the same way as it 
happens in “deformed QCD” model where all computa¬ 
tions are performed in a theoretically controllable regime 
[21-23]. 

In what follows we assume that the “strange” vacuum 
energy (6) and associated with this energy the topolog¬ 
ical susceptibility (10) is indeed saturated by fraction¬ 
ally charged monopoles with Q = ±l/iV which are con¬ 
stituents of KvBLL caloron with nontrivial holonomy 
(11). With this assumption in hand the question which 
is addressed in this work is as follows. How does the 
“strange energy” vary in a hyperbolic space x S;^_i 
as a function of dimensional parameter k? The differ¬ 
ence between the original x and x spaces 
is the curvature of the hyperbolic space i?[]Hlfi ^ at 
AC —>■ 0. If we find a linear dependence on ac at small ac 


it would be a strong argument supporting our conjecture 
(2) on linear dependence of “strange” vacuum energy as a 
function of external parameter. Such linear scaling obvi¬ 
ously implies that this background-dependent correction 
is not generated by any local operators such as curva¬ 
ture (4), but rather is generated by nonlocal operator 
(5), (11) which is sensitive to the global characteristics 
of the background. 


C. Nontrivial holonomy (11) in Euclidean space. 
Few technical details. 


In this section we present few formulae derived in Eu¬ 
clidean space in order to compare them with parallel ex¬ 
pressions obtained in the hyperbolic space x S],_i from 
section HI. The corresponding comparison will allow us 
to study the dependence of the “strange” energy as a 
function of At at small At —A 0. 

We start from analysis of the KvBLL configurations 
with nontrivial holonomy (11). We use SU{2) gauge 
group in our discussions to simplify notations, though 
the generalization for SU{N) is also known. The KvBLL 
caloron can be represented as a combination of two 
monopoles. The first monopole is a conventional BPS 
monopole, and at large separations between the con¬ 
stituents in Hedgehog gauge can be represented as fol¬ 
lows 



coth(ur) 


1 \ 

r) 2 r 



sinh(ur) 


iijk 


2r2 ’ 


( 12 ) 


where we adopted the notations from refs [28-34] to coin 
this constituent as M-monopole. Parameter v in this 
formula is arbitrary number which is determined by the 
holonomy (11). The classical moduli space is a circle, 

O-TT 

u G : 0 < u < —, (13) 


such that V is an angular variable. For any v ^ Q the 
gauge group is broken to t/(l). In other words, param¬ 
eter V plays the role of the vacuum expectation value of 
the Higgs field 4)(x), which is represented in this sys¬ 
tem by Al^-component of the gauge field as expressed 
by eq. (12). One should emphasize that we study glu- 
odynamics, without a scalar Higgs field in the system. 
Nevertheless, A 4 - component of the gauge field with non 
vanishing expectation value v plays exactly the same role 
as the Higgs field in adjoint representation in the stan¬ 
dard BPS equations. Therefore, it is not a mystery that 
the r dependence of the A 4 (ri)- component in eq. (12) is 
identically the same as in conventional BPS construction 
for the Higgs 4>(x) field®. 


® Our convention for normalization is: |<I>p = 2Tr('J>)'^ = 

2Tr [zl 4 (r —A 00 )]'^ = 






6 


The most important property of solution (12) is the 
behaviour of at large distances, which is convenient 
to represent in the unitary gauge: 

Af (r ^oo) ^ y, (14) 


where we also keep the Coulomb like correction 1/r along 
with the leading term ^ v because this Coulomb long 
range interaction will play an important role in our future 
discussions. One can explicitly see that the holonomy 
(11) is 


-TrL = cos 
2 



cos(7r:/), 


(15) 


where we introduced the dimensionless parameter v = 
1^. The holonomy belongs to the group center ^TrL = 
±1 when V assumes its boundary points {v = 0, ^), in 
which case it is called the trivial holonomy. 

The second type of monopole is the so-called L 
monopole which can be constructed from (12) as follows 
[26, 27]. First, one should replace u —^ — u, which is 
equivalent to replacement v ^ (1 — u). Than, one should 
make a “large” (improper) gauge transformation 

Ciarge(a;4) = exp ■ (16) 

As it is known the “large” gauge transformation should 
be treated differently from “small” (proper) gauge trans¬ 
formations because any two field configurations related 
by “large” gauge transformation do not belong to the 
same gauge orbit. Nevertheless, the transformation 
(16) preserves the periodic boundary conditions because 
fliarge(a ^4 = 0) = —[/iarge(a ^4 = P)■ Final Step is to per¬ 
form the reflection v —>■ —v in order to restore the original 
vacuum expectation value (13). It is implemented by the 
discrete transformation [/reflection = exp(fr^7r/2). The 
resulting configuration is the L monopole (the Kaluza- 
Klein monopole in the original terminology). Its asymp¬ 
totic behaviour is 

A^(r^oo)^ y, (17) 


which should be contrasted with (14) with an opposite 
sign for a Coulomb term. It corresponds to the opposite 
magnetic charges of the M and L monopoles. Therefore, 
the action S, topological charge Q, and magnetic charge 
q for M and L monopoles are: 

SM = ^v, Q^ = v, q^ = +l, i/=y, (18) 

5 27r 

= ^{1- ly), Q^ = {l-v), q^ = -l, 

r 

while the monopole’s mass m is determined as S' = m/3 
such that m = when m is expressed in terms of v. 
One can explicitly see from (18) that the classical action 
S = {S^ + S^) for the KvBLL configuration consisting 


L and M monopoles does not depend on v and coin¬ 
cides with action of the conventional periodic instanton 
[24, 25] with = 1 and action S = j. On 

the quantum level the partition function, of course, will 
depend on v. We review the relevant results from refs 
[28-34] on this matter below in section IID. Specifically, 
we want to pinpoint few crucial elements which differ be¬ 
tween Euclidean expressions and corresponding formulae 
written in hyperbolic space. Precisely this difference as 
we shall argue in section III is responsible for the linear 
in K correction in expression (2) describing the “strange” 
energy. 


D. The grand canonical ensemble of monopoles 
with nontrivial holonomy (11) 


In the semiclassical approximation the partition func¬ 
tion of the gluodynamics is represented by the statis¬ 
tical ensemble of an arbitrary number of interacting 
monopoles and anti-monopoles of all kinds. The cor¬ 
responding picture is well tested in the weakly coupled 
gauge theory, the so-called “deformed QCD” where all 
important elements such as the generation of the strange 
energy, the topological susceptibility, the contact term, 
etc, have been explicitly computed [21-23]. The key les¬ 
son from that studies is the crucial role of the holonomy 
(11) and non-locality in generation of all these effects. 
In fact one can argue [23] that the system belongs to a 
topologically ordered phase as a result of these non-local 
effects®. 

Similar computations in strongly coupled regime have 
not been performed yet. Nevertheless, one should expect 
a very similar behaviour of the strongly coupled ensemble 
of KvBLL calorons represented by a set of their monopole 
constituents as discussed above in section IIC. In other 
words, we assume that KvBLL calorons with nontrivial 
holonomy are responsible for generation of the strange 
energy, the topological susceptibility, the contact term, 
and many other highly nontrivial features of the system. 
One should emphasize that the corresponding contribu¬ 
tions are finite in the infrared (IR) in the large volume 
limit and small temperature, (V,P) —>■ oo, in contrast 
with conventional instanton computations. We do not 
claim that the semiclassical approximation adopted here 
is justified in the strongly coupled regime. In fact, it is ex¬ 
pected that the fluctuations with typical scales ^ Aqct 
change some numerical estimates. However, as we shall 
argue below the corresponding fluctuations with scales 


Unfortunately we can not use the “deformed QCD” model to 
address the question formulated in the present work. This is be¬ 
cause we can not implement k parameter into this model because 
the size of in “deformed QCD” model [21] must be small to 
keep the system in the weakly coupled regime, while in hyper¬ 
bolic space the radius of must be large as it is correlated 

with our large 3d volume 
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^ ^ QCD modify the contributions from the far 

IR regions with typical scales ^ k, which is precisely the 
subject of the present work. 

Therefore, we follow [28-34] and represent the grand 
partition function for the ensemble as follows 

Kl Km Ki Km 

E n n n n 

KlKmKj^Km iL = l iM = i iL = ^ = 1 

f fd?Xi^ f<PyrL fd-^y^M 

J Kl\ Km\ K-^\ K^\ 

X ■ detG[a;] • detG[y], (19) 

where / is the “fugacity” of the monopoles. The fugacity 
has been computed in [28] in terms of the fundamental 
parameters of the system such as A q^^ , /3, and holonomy 
V defined by eq. (13). One should emphasize that fugac¬ 
ity / is dimensional parameter which is sensitive to all 
scales of the problem as it is expressed in terms of zero 
as well as non-zero modes, see below. The 3 dimensional 
coordinates Xi ^, and describe the positions 

of M, L monopoles and L, M anti-monopoles correspond¬ 
ingly. The G[x\ is a (Kl + Km) x (Kl + Km) and G[y] is 
a {Ki + K^) X {Ki -I- Kj^) matrices describing the mod¬ 
uli space. Their explicit form are given in refs. [28, 29]. 
These matrices represent the standard zero mode con¬ 
tributions, and highly sensitive to the IR physics as they 
depend on the holonomy and the long range Coulomb in¬ 
teractions between the monopoles. The corresponding in¬ 
teractions can be traced from the asymptotical behaviour 
of the monopole’s solutions (14) and (17). Finally, poten¬ 
tial V{x—y) describes the interaction between monopoles 
and anti-monopoles of the entire ensemble. The corre¬ 
sponding interactions, along with G[x\ and G[y] are also 
long ranged. These elements of the partition function are 
highly sensitive to the IR behaviour and to the boundary 
conditions. 

There are many subtle points in writing (19) which 
shall not be discussed here. We refer to the original works 
for the discussions and references. For the moment we 
ignore the interaction between the monopoles and anti¬ 
monopoles, V{x — y). As we argue below in section III C, 
the corresponding interaction may change the numerical 
results, but can not modify our main claim on structure 
of the correction ~ k. 

The only relevant element for our future studies is the 
presence of the long range forces entering G[x\,G[y]. The 
corresponding interactions are effectively cancel in com¬ 
putation of the free energy due to the total neutrality 
condition^, as argued in [28, 29]. As a result of neu¬ 
trality the number of different types of monopoles is the 


^ The consequences on neutrality condition has been slightly cor¬ 
rected recently in [32—34]. The correct statement is not that 
the neutrality condition implies that the total charge is zero for 
each given configuration. Rather, the correct statement is that 
the corresponding expectation value of the charge vanishes while 
charge number density itself may still fluctuate. 


same in each given configuration. In other words, the 
partition function under these assumptions decouples for 
monopoles and anti-monopoles, Z = ■ Z-, and takes 

the following simple form [28, 29]: 


E 

KlKm 


(47r/F)^^+-^“ 

Kl\Km\ 




( 20 ) 


where V is the 3-volume of the system. The combina¬ 
tion v^^{l — which enters the partition function 

(20) comes form zero mode determinant which itself is 
expressed in terms of the classical actions of the con¬ 
stituents (18). In the large volume limit the sum is sat¬ 
urated by very large K such that the partition function 
can be evaluated using saddle point approximation. 

To proceed with estimations we first represent fc! using 
the Stirling formula 


1 


= e 




1 


-K\nK+K 


K\ y/2^ 

The next step is to replace the sum by the integrals 


( 21 ) 


^ dK dQ, 

KlKm '' 

K = Kl + Km, Q = Km-Kl, (22) 


where K describes the total number of monopoles in 
a given configuration, while Q describes the magnetic 
charge of a given configuration as M and L monopoles 
have opposite charges according to (18). Using saddle 
point approximation one arrives to the following expres¬ 
sion for the partition function in terms of the saddle value 
Kq which saturates (20) 


Z± ~ • J Ko = 8^/UvEi^-(23) 


Therefore the final expression for the partition function 
and the free energy for the vacuum ground state in this 
approximation assumes the form 

47rAl^ 

Z± ~ exp [47r/U], / = 

327r^ 

F_ = -rinZ = -^A|^U, (24) 

where we substitute the “confining” value for the holon¬ 
omy ly = 1/2 which minimizes the free energy. 

Few comments are in order. The expectation value 
{Q) = 0 obviously vanishes, such that system is neutral. 
However, the fluctuations of the do not vanish, in 
agreement with [32, 33], but strongly suppressed for large 
volume system, as expected 




1 

Vv 


0 . 


(25) 


The free energy (24) is finite at zero temperature limit. 
Furthermore, InZ is proportional to the 4-volume ^ V/T 


















of the system demonstrating the expected extensive scal¬ 
ing at low temperature. One should emphasize that these 
computations (based on configurations with nontrivial 
holonomy) generate the IR finite and well defined con¬ 
tributions to different observables expressed in terms of 
fundamental parameters of the theory, in contrast, for 
example, with instanton computations. The dynamically 
generated “confining” value for the holonomy zz = 1/2 is 
also a highly nontrivial phenomenon - it leads to a proper 
behaviour for the Polyakov’s loop, the Wilson loop and 
the string tension [28, 29], see also Appendix A with few 
historical and terminological comments on fractionally 
charged constituents, previously emerged in the litera¬ 
ture in different contexts and different systems. 

Furthermore, this “confining” value for the holonomy 
z/ = 1/2 leads to a consistent resolution of the so-called 
U{1)a problem formulated in terms of the topological 
susceptibility (8) and the 6 dependence of the “strange” 
vacuum energy (6). Indeed, the introduction of the 6 
term into the Lagrangian changes the fugacity for the 
monopole / —>■ and anti-monopoles / —>■ 

This modification follows from the fact that the topolog¬ 
ical charges for monopoles and anti-monopoles assume 
the magnitude Q = ±1/2 for “confining” value of the 
holonomy z^ = 1/2 as it follows from quantum numbers 
for monopoles (18). The anti-monopoles assume the op¬ 
posite sign for the topological charge as they are anti¬ 
selfdual solutions. This modification leads to replace¬ 
ment of expression (24) by the following formula which 
is valid for |0| < tt: 

Fva.c{S) = (^ 2 ) ■ 


The topological susceptibility now can be easily com¬ 
puted by differentiating (26) twice with respect to 6 with 
result 


X = 


1 d'^F^^ciO) _ a 4 


(27) 


Finally the vacuum energy (26) per unit volume 
FVac/I" ~ A4 is precisely the first term entering 
the expression (2). It has all the features of the 
“strange energy” briefly described in section IIA in 
model-independent generic way. The mechanism based 
on the KvBLL configurations reviewed above precisely 
generates all these required properties. Similar formulae 
can be easily generalized for arbitrary number of colours 
N when Fvac(^) ~ cos and x ^ 1, which is con¬ 
sistent with conventional resolution of the U(l) a problem 
in large N limit. 

We do not claim to have derived any new results in 
the present subsection. Rather, we just reproduced and 
explained the known results [28, 29, 32, 33] in slightly dif¬ 
ferent and simplified manner in order to analyze the role 
of similar vacuum configurations in cosmological context 
in next section III. Furthermore, we do not claim that the 
corresponding computations in strongly coupled regime 


are exact. In fact, we expect the corrections to be order of 
one to the fugacity / and all other numerical coefficients 
such as (24), (27) discussed above. However, we do not 
anticipate any drastic qualitative changes of this frame¬ 
work as a result of these possible corrections. In partic¬ 
ular, we expect that the free energy generated by these 
configurations remains finite in the IR and demonstrates 
the extensive behaviour at low temperature T —>■ 0 as 
presented above. Precisely these features will play a cru¬ 
cial role in our arguments on small modification of this 
“strange” vacuum energy with tiny variation of the back¬ 
ground to be considered in section HID. 


E. Interpretation. Cosmological context. 

There are many important elements related to the 
computations reviewed in previous sections. In what fol¬ 
lows we would like to make only very few comments which 
will be relevant for our studies in cosmological context on 
IR sensitivity of the system. 

1. First of all, the positive sign in (27) unambigu¬ 
ously implies that the corresponding configurations sat¬ 
urating the topological susceptibility (and related the 9- 
dependent portion of the vacuum energy) can not be 
identified with any propagating degrees of freedom in ac¬ 
cordance with (7). Indeed, the computations reviewed 
above explicitly show that the relevant configurations are 
the KvBLL calorons with nontrivial holonomy describing 
the tunnelling events between topologically distinct but 
physically identical winding states, rather than propa¬ 
gating gluons. All effects are obviously non-analytical in 
coupling constant ^ exp(—1/g^) and can not be seen in 
perturbation theory. 

2. One can view the relevant topological configura¬ 
tions as the 3d magnetic monopoles wrapping around 
time direction. This leads to the non-vanishing holonomy 
(11) and non-vanishing topological charge (18) of the con¬ 
stituents defined in 4d space-time. In the limit of T —>■ 0 
the “confining” value Z4 = 1/2 for the holonomy implies 
that the parameter v which determines the monopole’s 
mass in conventional 3d theory tends to zero in this limit, 
V = ttT —^ 0. At the same time the monopole gener¬ 
ates a finite contribution (24) to the path integral due 
to its finite 4d action (18) resulting from very long path 
~ T~^ in time direction. It would be misleading to inter¬ 
pret the confinement and other features (discussed above 
in section IID) generated by these configurations as the 
condensation of the monopoles. It would be more ap¬ 
propriate to use term “percolation” as the configurations 
described above do correlate at arbitrary large distances, 
but they obviously do not form a “condensate” in con¬ 
ventional condensed matter terminology. 

3. In the cosmological context such configurations are 
highly unusual objects: they obviously describe the non¬ 
local physics as the holonomy (11) is a nonlocal object. 
Indeed, the holonomy defines the dynamics along the en¬ 
tire history of evolution of the system in the given con- 
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fined phase: from the very beginning to the very end. 
There is no contradictions with causality in the system as 
there is no any physical degrees of freedom to propagate 
along this path at /3 —>■ oo, see item 1. above. Indeed, 
this entire gauge configuration is a mere saddle point 
in Euclidean (imaginary time) path integral computa¬ 
tion which describes the instantaneous tunnelling event, 
rather than propagation of a physical degree of freedom 
capable to carry an information/signal. 

4. Further to this point, the extensive property of free 
energy /3Fvac ^ E/T at /3 —>■ oo is a highly nontriv¬ 
ial phenomenon in this framework as the 4-volume VjT 
appears in this description due to few important steps. 
First, one should regularize the moduli space by cutting 
off an each given configuration in the IR. Secondly, one 
should sum over all configurations (20) by using saddle 
point approximation in large volume limit, which even¬ 
tually leads to (24). This “emergent” extensive property 
is drastically different from conventional approaches to 
cosmology when the free energy is determined by the 
Lagrangian density T[4>] integrated over the 4-volume 
J d^x. In this last case the extensive property is a triv¬ 
ial manifestation of the system formulated from the very 
beginning in terms of the local field 4>(a;). It shows one 
more time that generation of the “strange” energy (2) 
is highly non-local non-perturbative effect when the vol¬ 
ume of the system could be very large, but still finite to 
proceed with computations (24) in this approach. Essen¬ 
tially, the finite local energy density of the system (24) in 
this framework is determined by the entire time evolution 
/3 —>■ oo in confined phase, which is obviously a non-local 
procedure. Still, it does not contradict the causality, see 
item 3 above. 

5. Last, but not least. All these highly nontrivial 
non-local features listed above emerge only at T < Tc 
when the configurations with nontrivial holonomy (11) 
start to play the dominant role in the dynamics. At high 
temperatures the contribution of the configurations with 
nontrivial holonomy can be completely ignored as they do 
not contribute to the partition function in thermodynam¬ 
ical limit. This property of drastic variation of “strange” 
energy (2) with temperature around Tc may play an im¬ 
portant role in cosmological context as we discuss in next 
section. 


III. NON-TRIVIAL HOLONOMY AND 
HYPERBOLIC SPACE x 

The main goal of this section is to generalize the re¬ 
sults presented above in sections IIC, IID to hyperbolic 
space x to argue that the correction to the free 
energy (24) are linearly proportional to k at small k —>■ 0. 
In this limit it is quite obvious that all features listed in 
section IIE on nature of the “strange energy”, including 
its non-local nature, remain the same in this limit k —>■ 0 
as the system is almost 4d Euclidean space, with very 
tiny deviations ~ k which we wish to recover. The corre¬ 


sponding linear dependence on k would strongly support 
our conjecture that the correction in eq. (2) are linearly 
proportional to the Hubble constant®, which dynamically 
drugs our Universe to the de Sitter state as Friedman 
equation (3) suggests. 

We discuss the relevant gauge configurations in hyper¬ 
bolic space in subsection III A, while the grand canonical 
ensemble of such hyperbolic monopoles will be studied in 
subsection III B where we discuss the crucial distinctions 
between hyperbolic and Euclidean monopoles. Precisely 
this difference eventually leads to a tiny ^ k variation 
from the Euclidean results. The corresponding devia¬ 
tion is expressed in terms of the fugacity in section III C, 
which ultimately leads to slight modification of the vac¬ 
uum energy. We list some profound cosmological conse¬ 
quences of this modification of the vacuum energy with 
background in section HID. 


A. Holonomy and monopoles in hyperbolic space 


The construction of the monopoles on hyperbolic space 
has been known since [35-37]. Furthermore, many 
other topological objects, including calorons, instantons, 
vortices, skyrmions have been constructed on hyperbolic 
space [38-43] . An important technical element which was 
used in these constructions is the conformal equivalence 
of and x S].-i. Indeed, this equivalence can be 
explicitly checked by introduction toroidal coordinates 
{p, 6, (j), x) on as follows: 


Xi 


X2 


x^ 


Xi 


(a:i,a;2,a:3,X4) 


I 

cosh(Kp) -I- cos X 

I 

cosh(Kp) -I- cos X 

I 

cosh(Kp) -I- cos X 


sinh(Kp) sin 9 cos </>, 
sinh(Kp) sin 0 sin (j), 
sinh(Kp) cos 9, 


I 


cosh(Kp) -|- cos X 


sinx. 


(28) 


It is then easy to check that the metric on R^ becomes 

ds^(]HI®) -I- K~‘^dx^ 


ds^ 


)= dx^dx^ = 


(cosh(«:p) -|- cos x)^ ’ 


(29) 


where (is^(]HI®) is the metric on hyperbolic 3-space with 
spherical coordinates (p, 0, <()): 


ds^(]HI®) = dp^ + 0d(j)^) . (30) 

The holonomy in terms of these variables is computed 
along the circle §].-i parameterized by X; that is 

U = Vexp(^^J dxA^^ , (31) 


see footnote 3 with some clarification on terminology. 
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where is the component of the gauge potential associ¬ 
ated with coordinate x, and it plays the same role as the 
A 4 , similar to the construction in Euclidean space (12). 
In both cases, the A 4 and A^ assume non-vanishing ex¬ 
pectation values, and play the same role as the Higgs field 
<i> in adjoint representation, as explained after eq. (13), 
see also footnote 5 on our normalization. Formula (31) 
plays the same role as equations (5), (11), while dx/ k, in 
eq. (31) plays the role of dx 4 in evaluation of the holon- 
omy computed along the circle according to (13). At 
large p —>■ 00 the approaches a non-vanishing constant 
value, similar to parameter v in Euclidean space (12). 

With these remarks in mind the explicit form for the 
Bogomolny- Prasad- Sommerfeld (BPS) monopole in hy¬ 
perbolic space in the unitary gauge can be written as 
follows [36, 37]: 

t-3 

A^ (p) = [C'KCoth(C'Kp) — KCoth(Kp)] (32) 

where parameter C takes any value greater than 1. In 
formula (32) we limited ourselves by writing down only 
the A^^(p)-component, which defines the boundary con¬ 
ditions at large p. We coin this solution as M monopoles 
in order to be consistent with the terminology introduced 
for the Euclidean counterparts (12). 

While many topological objects, including calorons 
with trivial holonomy [41] have been constructed in hy¬ 
perbolic space, as we already mentioned, an explicit con¬ 
struction of the calorons with nontrivial holonomy, which 
is analogous to the KvBLL solutions [26, 27], has not 
been constructed yet. In what follows we shall assume 
that such solutions do exist, though we do not need their 
explicit form in our future discussions. Important point is 
that if such configurations exist than they must exhibit 
the same features which the KvBLL solutions demon¬ 
strate. Namely they could be viewed as a set of 2 differ¬ 
ent types of monopoles for SU{2) gauge group. The first 
type is precisely the M-monopole (32) discussed above, 
while the second one, the L monopole can be constructed 
as described in section IIC for the Euclidean counterpart. 
We shall return to this construction later in the text, but 
first, we want to understand the physical meaning of the 
parameter C entering eq. (32) by analyzing the limit 
At —>■ 0 when the Euclidean monopole (12) is recovered. 

To recover the Euclidean monopole solution one should 
take the limit C —> 00 along with /t —> 0 with combina¬ 
tion Ck being fixed to be equal v. In this limit A^ (p) 
becomes 


for A^ (p) and holonomy 

[^(C ~ 1) + (^(e ^)] -^1 

iTr{7(p — 00 ) = cosTTi/, i/=(C'—1), (34) 

where we introduced parameter iz expressed in terms of 
original parameter C. It plays the same role as parameter 
iz discussed in the Euclidean construction (5), (18). The 
crucial observation here is that the asymptotic formula 
for A^ {p) does not contain a long range Coulomb in¬ 
teraction p~^ in contrast with its Euclidean counterpart 
(14). Instead, there is an exponentially suppressed cor¬ 
rection ~ exp(—/tp) in formula (34). One can interpret 
such drastic changes in behaviour of the solution as a re¬ 
sult of screening of the magnetic field by the curvature 
in the hyperbolic space. In terms of the new parameter 
iz the solution (32) assumes the form 

3 

A^{p) = ^(zz -t-1) coth [(zz -I- 1)kp] — cothrep^ ^^.(35) 

Our next step is to recover the L monopole assuming 
that KvBLL caloron with nontrivial holonomy in the hy¬ 
perbolic space exists, similar to construction [26, 27] in 
the Euclidean space. We follow the same steps to re¬ 
construct their properties presented in section IIC for 
the Euclidean monopolies. The first step is to replace 
zz —>• (I — zz). The second step is to make a “large” gauge 
transformation which assumes the following form in hy¬ 
perbolic variables 

Hiarge(x) = oxp (*yx) ■ (36) 

As we already mentioned the “large” gauge transforma¬ 
tions should be treated differently from “small” (proper) 
gauge transformations because any two field configura¬ 
tions related by “large” gauge transformation do not be¬ 
long to the same gauge orbit. Nevertheless, the trans¬ 
formation (36) preserves the periodic boundary condi¬ 
tions for the fields in the adjoint representation because 
C^iarge(x = 0) = -C/iarge(x = 27r). The final step 
is to perform the discrete transformation I/reflection = 
exp(iT^7r/2) to restore the original boundary conditions 
(34). The resulting configuration is the hyperbolic L- 
monopole: 

3 

A^{p) = ^1 — (zz -I- 1) coth [(zz -b 1 )kp] + cothztp^ ^^,(37) 


K^p) = 


V coth(up) — 


Ck = V (33) 


where p should be identified with r in the Euclidean 
space. Expression (33) exactly coincides with (12) and 
its asymptotic behaviour (14) in the unitary gauge. 

We are now in position to discuss the asymptotical 
behaviour of (32) and the holonomy (31) for finite k. 
Taking p —>■ oo one arrives to the following expressions 


where we introduced zz = (1 — zz) for convenience. One 
should emphasize that electric and magnetic fields of the 
the L -monopoles do depend on x variable as a result 
of X dependent “large” gauge transformation (36). This 
is analogous to L -monopole solution in Euclidean space 
(17) which is a time dependent configuration, and cease 
to exist in static 3d space. 

The classical action, topological and magnetic charges 
of the M and L constituents are determined by the 
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boundary conditions (35) and (37) at large p —>■ oo, sim¬ 
ilar to the Euclidean counterparts (18). The correspond¬ 
ing parameters S,Q,q obviously assume the same values 
(18) when expressed in terms of ly. These dimensionless 
parameters obviously can not depend on dimensional pa¬ 
rameter K, including k —>■ 0 limit. In fact, the correspond¬ 
ing formula relating Q and q for BPS M- type monopole 
(35), which identically coincides with the Euclidean ex¬ 
pression, was derived in hyperbolic space in the original 
work [35] , while the relation between S and Q is a direct 
consequence of self -duality of Yang-Mills equations. The 
monopole’s mass (the total energy of the configuration), 
being a dimensional parameter, does depend on k. Mass 
m satisfies an obvious relation S = m^, similar to its 
Euclidean counterpart, see text after eq. (18). Explic¬ 
itly, in terms of the original parameter C it is given by 
m = ^k{C — 1), which reduces to its Euclidean form 
m = when one takes the corresponding limit (33). 


B. Grand partition function for monopoles in the 
hyperbolic space 

With our main assumption that the calorons with non¬ 
trivial holonomy exist in hyperbolic space, similar to Eu¬ 
clidean KvBLL construction [26, 27], one should expect 
that the corresponding grand partition function Z has 
the following form, which is analogous to expression (19) 
discussed in previous section; 

Kl Km Kj^ Km 

E n n n n 

KlKmK[^Km ii, = liM = li£ = l*M = l 

^ f fy/gtPxiM lygd^ViL 

As Kl'. Km'. K^\ Km'. 

X ■ detG[a;] • det G[y], (38) 

where / is the “fugacity” of the hyperbolic monopoles. 
The corresponding dimensional parameter is highly sen¬ 
sitive to many details of monopole’s structure and their 
interactions with other monopoles. It is obviously differ¬ 
ent from its cousin fugacity / computed in the Euclidean 
space as discussed in section IID. We shall estimate the 
corrections to / later in section III C where we argue that 
the difference (/ — /) ^ k is linear in k at small k —>■ 0 
for configurations with nontrivial holonomy. 

The next item to discuss from formula (38) is the 3— 
dimensional coordinates XiM, and , yi^. They de¬ 
scribe the positions of M, L monopoles and L, M anti¬ 
monopoles correspondingly. These coordinates play the 
same role as in formula (19) with the only difference is 
that the distance between constituents is computed using 
the metric (30). The G[a;] is a (AY -|-ATm) x (ATl -I-ATm) 
and G[y\ is a {Ki + Km) x [Ki + Km) matrices de¬ 
scribing the moduli space. The difference with corre¬ 
sponding Euclidean expressions is that the behaviour at 
large distances is not the Coulomb like, but rather the 
exponentially suppressed as we already mentioned (34). 


This is because the matrices G[a;] and G[y] are computed 
from the corresponding zero modes in the background of 
the monopoles. At the same time, the zero modes, as 
usual, are fixed by the corresponding classical solutions. 
Therefore, the asymptotical behaviour of classical solu¬ 
tions (35) and (37) dictates the behaviour of G[x\ and 
G[y\. In both cases the computations of the free energy 
is reduced to expression (20) as a result of neutrality con¬ 
dition discussed in section IID. 

Final item to discuss is the factor ^ which accounts 
for the curvature of the hyperbolic space H^. It enters 
along with the spatial volume of the system where com¬ 
putations are being performed. We put the system into 
a large volume of radius R such that the volume is 

^ sin 0 sinh^ Acp 
nRy V5=- - -^.(39) 


47r f sinh2KA 

^ V 2 


It reduces to the Euclidean expression V = in the 
limit of small k with the corrections of order 0{k^) which 
are consistently neglected in the present work. Therefore, 
by repeating all the steps leading to formula (24) we ar¬ 
rive to the following expression for the partition function 
and the free energy in hyperbolic space 

Z± ~ exp [47r/E] , Avac = InZ, (40) 


where we identify from x geometry (reviewed 
in section II) with circle from x geometry 
presented in section III A. In other words, we identify 



K 


(41) 


such that the corresponding expressions for holonomy 
(5) and (31) coincide. With this identification the cor¬ 
responding formulae (15) and (34) in terms of dimen¬ 
sionless parameter i/ also coincide. In formula (40) we 
substitute the confining value for holonomy i/ = 1/2 as 
it has been done in the Euclidean space. This is because 
the free energy is minimized at zz = 1/2 irrespectively to 
the value of the fugacity, which is indeed different for two 
different geometries. To conclude: the only difference be¬ 
tween Z and Z describing the system on x S/_i and 
X geometries correspondingly is that the fugacities 
in these two systems assume slightly different values at 
small K —?► 0, which is the subject of the next subsection. 


C. Monopole’s fugacity 

We start our analysis with explanations on how the 
basic dimensional parameter, the fugacity, emerges in 
the system. This x independent dimensional parame¬ 
ter / effectively determines the dynamics of the system. 
This parameter essentially represents the density of the 
monopoles in the system. The classical action, the zero 
and nonzero mode contributions lead to the following ex¬ 
pression for the monopole’s fugacity / in terms of the 
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fundamental parameters of the theory [28, 29]: 


47r/3Ai 


QCD 


f = 


[1 + 2Trvu^] 

= = (77-^li 

(^QCD ^12) 




(42) 




where brackets (...) imply averaging over separation ri 2 
between M and L monopoles in ensemble (19) such as / 
is a;-independent as it should. 

Few comments are in order. Each KvBLL caloron is 
represented by the L and M monopoles and accompanied 
by 8 zero modes. It explains the major dimensional factor 
in eq. (42), including . The remaining numerical 

dimensionless factor “c” entering (42) is order of 1. It 
includes factor (A q^j-, )"^/^ which can be easily restored 
from the renormalization group analysis which requires 
that Aq^jj enters (42) with power (Aq^)^^/^. Subse¬ 
quently, this factor (A q^q )"^/^ must be accompanied by 
a dimensional parameter, which at small temperatures 
could be nothing else but the separation distance 
between the monopoles. In estimates of refs. [29] the nu¬ 
merical coefficient c is assumed to be one, which precisely 
corresponds to the expression for / given in eq. (24). 

Next factor which can be easily explained is the first 
term in numerator, [I -|- 2'kvv^]. This term has been 
originally computed in [26] and reproduced in [28, 29]. It 
is originated from the zero mode determinant. The cru¬ 
cial point here is that the algebraic dependence on ri 2 
emerges as a result of long range Coulomb terms in the 
classical solutions (14) and (17). This is because the zero 
mode structure is unambiguously fixed by the classical so¬ 
lutions with the corresponding Coulomb terms. Another 
important element is that this terms is proportional to 
the holonomy vv. It implies that this term will not be 
generated for the configurations with trivial holonomy. 

The nature of next two terms in numerator in eq. (42) 
is much harder to explain because they are originated 
from the contributions of the nonzero modes. The only 
comment we would like to make here is that these terms 
are also proportional to the holonomy, and can not be 
generated by the configurations with trivial holonomy. 

Now we are in position to estimate the difference be¬ 
tween the fugacity generated by monopoles in the Eu¬ 
clidean space versus hyperbolic monopoles at small k —>■ 
0. The estimation is convenient to represent in the fol¬ 
lowing form 


^ ~ [1 + Azm) • (1 -|- Anzm), A^m ~ -) (43) 

J ^ ^''QCD 

where factors A^m and Anzm describe the corrections due 
to the zero and non-zero modes correspondingly. The 
correction factor Azm comes from the first term in nu¬ 
merator (42), where we estimate (ri 2 ) ^ as the 

only dimensional parameter in the system at small tem¬ 
perature. We also expressed the result in terms of k 


rather than /3 according to identihcation (41). The cru¬ 
cial point in our estimate Azm is that a similar correction 
in hyperbolic space is absent as the corresponding classi¬ 
cal solutions (35) and (37) have exponentially suppressed 
asymptotic at large distances, in contrast with the Eu¬ 
clidean counterpart. Therefore, a term ^ ri 2 can not 
be generated in hyperbolic space, in contrast with the 
Euclidean case (42). 

Unfortunately, a similar unambiguous conclusion can 
not be reached regarding the nonzero mode contribution 
Anzm- This is because both, the Euclidean as well as hy¬ 
perbolic monopoles may generate such contributions pro¬ 
portional to K ■ ri 2 , which eventually produce a desired 
correction ^ k/Aq^. While the Euclidean expression 
is known and is represented by two factors in numera¬ 
tor in eq. (42), a similar expression in hyperbolic space 
is not known simply because an explicit construction of 
the hyperbolic calorons with nontrivial holonomy, similar 
to KvBLL solution, is yet unknown. Therefore, for the 
numerical estimates in what follows we set Anzm = 0- 

Few comment are in order. First, we want to argue 
that an unknown Anzm correction can not exactly cancel 
the computed Azm term (43). Indeed, the Azm structure 
is precisely fixed by the structure of the SU{N) gauge 
group with 4iV zero modes, while Anzm varies and de¬ 
pends, in particular, on presence of the matter helds, 
and other details of the system. In other words, a pos¬ 
sible cancellation, if ever occurs, can not be a universal 
phenomenon. Therefore, we use Azm and disregard Anzm 
as our order of magnitude estimate for ratio (43) . It ex¬ 
plicitly exhibits the linear in k correction to the fugacity, 
which is the main result of the present work. This cor¬ 
rection can be only generated by the configurations with 
nontrivial holonomy. In particular, conventional instan- 
tons and calorons with trivial holonomy may only gen¬ 
erate the higher order corrections ^ and do not 

contribute to the linear term (43). 

We should emphasize that estimate (43) was derived 
under assumption that the interaction term V{x — y) 
in formula (38) vanishes. It corresponds to the en¬ 
semble containing exclusively the monopoles (or anti¬ 
monopoles). Only in this case the partition function is 
exactly reduced to simple form (20) as a result of neu¬ 
trality condition as argued in [28, 29]. In reality the in¬ 
teraction plays crucial numerical role at finite temper¬ 
ature as shown in [32-34]. However, our main claim 
is that the linear correction (43) may receive large nu¬ 
merical corrections, but it can not be exactly cancelled 
as a result of unaccounted interaction term V{x — y) in 
formula (38). The basic argument behind this claim is 
the same one as presented above and based on observa¬ 
tion that the interaction V{x — y) is highly sensitive to 
the matter content of the theory (the number of flavours 
and its masses in the system), while a Azm in formula 
(43) is not sensitive to these modifications. The interac¬ 
tion V(x — y) may change a numerical coefficient in the 
estimate (ri 2 ) ^ A^^ which enters (43), but can not 
completely destroy this term. Therefore, our main claim 
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(that the linear correction (43) will be generated) holds 
irrespectively to any type of monopole-anti monopole in¬ 
teractions V(x — y). 


D. Linear corrections ~ k to the vacuum energy 


The result (43) for monopole’s fugacity can be trans¬ 
lated into the statement on variation of the vacuum en¬ 
ergy density in the bulk of space-time with a tiny varia¬ 
tion of the background. Indeed, according to (24), (40) 
and (43) the relevant ratio for the vacuum energies at 
K —>■ 0 for two different geometries can be represented as 
follows 

X J ( VV K 

i;vac[K3x§i] -J- [ ~Y’a= 

The same result can be represented in more conventional 
form 



SvacP^ X Sl_i] -A^ 1 


~ -A: 


QCD 


-I- K • A^ 


QCD 


QCD 2 ’ 


T 



(45) 


where we omitted all irrelevant numerical factors in ex¬ 
pression for the vacuum energy in Euclidean space, but 
kept the relevant sign minus (—) in front, which is well 
known feature of QCD. Our final formula (45) is a precise 
analog (in a simplified model) for the vacuum energy (2) 
conjectured for the de Sitter space. As we emphasized in 
the Introduction, the significance of the linear correction 
in eq. (2) is that the Friedman equation (3) unambigu¬ 
ously predicts a non trivial solution with constant Hq 
if the subtraction procedure is adopted as discussed in 
Introduction. The constant solution Hq automatically 
corresponds to a desired de Sitter behaviour (1), which 
might be relevant for the early Universe during the in¬ 
flationary epoch, and in present epoch for description of 
the dark energy. Few comments are in order: 

1. The difference between two geometries, x S),_i 
and x when the sizes of S^_i and are identi¬ 
cally coincide according to (41) is the small curvature 
^ of the hyperbolic space. According to conventional 
arguments on locality as discussed in Introduction it un¬ 
ambiguously suggests that all corrections must be pro¬ 
portional to the even powers However, we obviously 
observe a liner correction (45) in explicit computations. 

2. This linear correction ~ is generated by the config¬ 
urations with nontrivial holonomy, which itself is a non¬ 
local, but gauge invariant operator. Therefore, the stan¬ 
dard arguments on locality, reviewed in Introduction are 
badly violated by such configurations. One can see from 
eq. (45) that the linear correction ^ k is explicitly pro¬ 
portional to the holonomy i/, which is the gauge invari¬ 
ant observable, not reducible to the curvature. In other 
words, this correction is generated by non-local configu¬ 
rations, and can not be expressed in terms of local cur¬ 
vature ^ 


3. All effects discussed in the present work are non- 
analytical in coupling constant ~ exp(—I/g^) and can 
not be seen in perturbation theory. 

4. The result (45) is consistent with the previous anal¬ 
ysis in weakly coupled “deformed QCD” model where 
one can study the sensitivity of the vacuum energy to 
the very large distances by putting the system into the 
box of size L. It turns out [44] that the corrections to 
the vacuum energy are linear in inverse size ^ L“^. This 
model is very similar in all respects to the system studied 
in the present work because the vacuum energy in “de¬ 
formed QCD” model is also saturated by the monopoles 
with nontrivial holonomy. At the same time the con¬ 
ventional instantons with trivial holonomy produce only 
quadratic corrections ^ as noticed in [44]. 

5. The generation of the linear correction ~ k is also 
consistent with computations of ref. [4] in “deformed 
QCD” model where analysis was performed in terms of 
auxiliary topological non-propagating field. In that com¬ 
putations the root of the phenomenon is the presence 
of the non-trivial holonomy and long range monopole’s 
field which eventually is responsible for generation of the 
linear correction. 

6. The linear correction observed in our work is also 
consistent with the lattice simulations [45] when one 
studies the dependence of the vacuum energy on the size 
of the system. 

7. Our results are also consistent with the lattice sim¬ 
ulations [46] when the author studies the rate of particle 
production in the de Sitter background. The rate turns 
out to be linearly proportional to the Hubble constant 
~ H, rather than H^. Our comment here is that the 
rate of the particle production in quantum field theory 
in general is determined by the imaginary part of the 
stress tensor, Im[T)(], while the vacuum energy is related 
to the real part of the stress tensor, Re[T^]. Analytic- 
ity suggests that both components must have the same 
corrections on H at small H. Therefore, the lattice mea¬ 
surements [46] of the linear dependence on H strongly 
suggest that the vacuum energy (which is determined by 
the real part of the same stress tensor) must also ex¬ 
hibit the same linear ^ H correction. The corresponding 
lattice computations of the 0 dependent portion of the 
vacuum energy and topological susceptibility in time de¬ 
pendent background are possible in principle, but tech¬ 
nically much more involved than the analysis performed 
in ref. [46]. 

8. Last but not least. The sign for the difference 

AFvac = X S^_i] - E[gij,^ = is positive as 

one can see from eq. (45). It corresponds to the posi¬ 
tive sign for the cosmological constant (dark energy) in 
cosmological context. 

There is a fundamental difference in signs with conven¬ 
tional Casimir effect when the corresponding subtraction 
procedure typically leads to the negative, rather than 
positive, sign for the vacuum energy. This difference is 
due to the fact that the conventional Casimir vacuum 
energy is generated by the fluctuations of the physical 
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propagating photons. It is drastically different from the 
the vacuum energy computed in the present work when it 
is generated by the tunnelling transitions between differ¬ 
ent topological sectors. As we explained in section IIA 
the corresponding vacuum energy can not be expressed 
in terms of any propagating degrees of freedom as it has 
pure non-dispersive nature. This is precisely the origin 
for positive sign of the vacuum energy: AAvac > 0. 


IV. CONCLUSION 

The formal result of the present work can be expressed 
by eqs.(44), (45), and we shall not repeat the comments 
listed in the last section III D explaining some important 
consequences of this result. If the same effect persists in 
FLRW Universe (2), which we expect to be the case, it 
may have a number of profound consequences for under¬ 
standing of the past, present and future evolution of our 
Universe. 

First of all, the nontrivial holonomy (5) implies the 
presence of as a part of our space-time of our Uni¬ 
verse. It is an additional invariant characteristic of the 
manifold which can not be reduced to the local curva¬ 
ture. In construction discussed in the present work the 
corresponding is identified with Euclidean time direc¬ 
tion. This is because the original intention in early works 
on the subject (where the corresponding topological vac¬ 
uum configurations were invented) was to analyze the 
temperature dependence of the QCD phase transition. 
In principle, similar could be also a part of spatial co¬ 
ordinates. Such an assumption definitely consistent with 
all known observations if the size of the is sufficiently 
large > H~^ at present epoch, see [47] for the estimates 
in the given context. Furthermore, the linear correction 
enters formula (45) in form of its absolute value |«:| as it 
essentially describes the linear (positively defined) size of 
the corresponding manifold. In context of FLRW Uni¬ 
verse a similar statement implies that the linear correc¬ 
tion in H enters formula (2) in form of its absolute value 
\H\, see footnote 3 for terminological clarification. There¬ 
fore, it can not lead to any T-violating effects, which one 
could suspect as H = aja indeed is a T-odd parameter. 
Still, it generates the de Sitter behaviour (1) as a result 
of this linear scaling. 

Another profound consequence of this framework is as 
follows. The conventional scenarios of the eternal self- 
producing inflationary universes are always formulated in 
terms of a physical scalar dynamical inflaton field 4>(a;). 
This problem with self-reproduction of the universe does 
not even emerge in our framework as there are no any fun¬ 
damental scalar dynamical fields in the system responsi¬ 
ble for inflation. Instead, the de Sitter behaviour (I) in 
our framework is pure quantum phenomenon, which is 
a consequence of the dynamics of the long ranged topo¬ 
logical configurations with nontrivial holonomy, rather 
than a result of a physical fluctuating dynamical field. A 
“strange nature” of this type of energy manifests itself 


in terms of the “wrong” sign in the correlation function 
which can not be formulated in terms of any local prop¬ 
agating degrees of freedom as explained in section II A. 
The corresponding topological configurations which are 
responsible for this behaviour may generate, as argued 
in this work, the linear in H correction in the Friedman 
equation (3) which eventually leads to the de Sitter be¬ 
haviour. Other problems formulated in terms of scalar in- 
fiaton field <I’(a:) (such as large initial value <I>in Mpl 
for the inflaton) do not emerge in this framework, see 
[4, 5] for the details. 

Finally, we should mention that the energy described 
by a formula similar to eq. (2) (which eventually leads 
to the de Sitter behaviour (1)) has been previously pos¬ 
tulated [48, 49] as the driving force for the dark energy. 
The model has been (successfully) confronted with obser¬ 
vations, see recent review papers [50, 51] and many origi¬ 
nal references therein, where it has been claimed that this 
proposal is consistent with all presently available data, 
see also ref. [52] for completeness. Our comment here is 
that history of evolution of the universe may repeat itself 
by realizing the de Sitter behaviour twice in its history. 
The QCD-dynamics was responsible for the inflation in 
early universe, while the QCD dynamics is responsible 
for the dark energy in present epoch. 

We conclude this work (mainly devoted to analysis of 
the topological configurations with typical energy scale 
A qcd ) with the following comment related to a funda¬ 
mentally different problem with drastically different en¬ 
ergy scale. Namely, as we discussed at length in this 
paper, the heart of the proposal is a fundamentally new 
type of energy (2), (24), (45) which can not be expressed 
in terms of any propagating degrees of freedom. Rather, 
this novel contribution to the energy has non-dispersive 
nature. The effect is formulated in terms of the tunnelling 
processes between topologically different but physically 
identical states. This novel type of energy, in fact, has 
been well studied in the QCD lattice simulations in the 
flat background, see [4] for references on the original lat¬ 
tice results. Our comment relevant for the present study 
is that this fundamentally new type of energy can be, in 
principle, studied in a tabletop experiment by measur¬ 
ing some specific corrections to the Casimir vacuum en¬ 
ergy in the Maxwell theory as suggested in [53-56] . This 
fundamentally new contribution to the Casimir pressure 
emerges as a result of tunnelling processes, rather than 
due to the conventional fluctuations of the propagating 
photons with two physical transverse polarizations. This 
effect does not occur for the scalar field theory, in con¬ 
trast with conventional Casimir effect which is opera¬ 
tional for both: scalar as well as for Maxwell fields. The 
extra energy computed in [53-56] is the direct analog of 
the non-dispersive contribution to the energy (2) which 
is the key player of the present work. In fact, an extra 
contribution to the Casimir pressure emerges in this sys¬ 
tem as a result of nontrivial holonomy similar to (5) for 
the Maxwell field. The nontrivial holonomy is enforced 
by the nontrivial boundary conditions imposed in refs 
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[53-56]. 
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Appendix A: Few comments on fractionally charged 
constituents and the terminology 

The constituents of the KvBLL conhgurations were 
originally [26, 27] called the BPS (Bogomolnyi-Prasad- 
Sommerfeld) monopoles (14) and KK (Kaluza- Klein)- 
monopoles (17) correspondingly. These configurations 
were later on coined as M- dyons and L- dyons to empha¬ 
size that they carry the electric charges along with the 
magnetic charges [28, 29]. This (incorrect) interpreta¬ 
tion based on observation that these configurations carry 
the topological charges (and naively the electric charges) 
along with magnetic charges. One should remember, 
however, that the monopoles in this construction are 
pseudoparticles living in 4d Euclidean space-time, rather 
than static 3d objects. The finite action and finite topo¬ 
logical charge for these objects results from wrapping of 
the monopole’s path along the Euclidean time direction 
with nontrivial holonomy (11). These objects do not 
carry a conventional static electric charge; nevertheless, 
they do carry the topological charges defined in 4d Eu¬ 
clidean space-time. Furthermore, the second types of the 
monopoles, the L monopoles are time dependent configu¬ 
rations, and do not exist as static objects in 3d Euclidean 
space. Therefore, we keep notations for letters M and L 
suggested in [28, 29], but we use term “monopoles” rather 
than “dyons” in the present work. When the holonomy 
assumes its “confining” value v = \/2 the topological 
charges of the constituents assume Q = ±1/2. 

In more generic case of the SU{N) gauge group the 
topological charge Q = zLl/N for “confining” holonomy 
such that a single KvBLL configuration can be thought as 
a superposition of N different types of monopoles which 
carry N different types of magnetic charges and fractional 
1/N topological charge such that the superposition car¬ 
ries an integer topological charge. 

We believe a short historical detour on fractionaliza- 
tion of the topological charge in QFT is warranted here. 
In the given context fractional topological objects appear 
in 2 dimensional CP^~^ model [57] which were coined as 
instanton quarks (other names: instanton partons, frac¬ 


tional instantons). These quantum objects carry frac¬ 
tional topological charge Q = ±1/N, and they are very 
similar to the L and M monopoles discussed in this work. 
These objects do not appear individually in path integral; 
instead, they appear as configurations consisting N dif¬ 
ferent objects with fractional charge 1/A^ such that the 
total topological charge of each configuration is always 
integer. In this case 4Afc zero modes for k instanton so¬ 
lution is interpreted as 4 translation zero modes modes 
accompanied by every single instanton quark. The same 
counting holds, in fact, for any gauge group G, not lim¬ 
ited to SU{N) case. While the instanton quarks emerge 
in the path integral coherently, these objects are highly 
delocalized: they may emerge on opposite sides of the 
space time or be close to each other with alike probabil¬ 
ities. Similar objects have been discussed in a number 
of papers in different contexts, including the topic of the 
present work, [21, 22, 26-29, 58-63]. 

In particular, it has been argued that the well- 
established 9/N dependence in strongly coupled QCD 
(expressed by formula (26) for specific case N = 2) un¬ 
ambiguously implies that the relevant configurations in 
QCD must carry fractional topological charges in confine¬ 
ment phase, see review preprint [59] and the references 
on earlier original results therein. The weakly coupled 
deformed QCD model [21-23] where computations are 
under complete theoretical control is a precise dynami¬ 
cal realization of this idea when the fractionally charged 
monopoles are responsible for confinement, saturate the 
topological susceptibility with a “wrong sign”, generate 
the “secrete long range forces”, suspected long ago [20], 
and provide other crucial elements which are known to 
exist in strongly coupled regime as reviewed in section 
II A. 

Furthermore, it has been argued in [60, 63] that the 
confinement deconfinement phase transition within this 
framework can be interpreted as Berezinskii-Kosterlitz- 
Thouless (BKT) -like phase transition: at T > Tc the 
constituents prefer to organize a single caloron of a fi¬ 
nite size. We coin this phase as a “molecular phase” 
which corresponds to a de-confined phase in conven¬ 
tional terminology. When one crosses the phase tran¬ 
sition line at T < Tc the constituents (which are called 
L, M monopoles in the present work) prefer to stay far 
away from each other. It corresponds to the dissociation 
of each caloron into N constituents, and we call this state 
as the “A component plasma phase” in 4d Euclidean 
space. This regime corresponds to the confined phase in 
conventional terminology when all constituents are delo¬ 
calized in 4d Euclidean space. The gap in this confined 
phase is determined by the Debye correlation length of 
this 4d plasma. The arguments [60, 63] are based on 
large N counting, but we believe that this picture holds 
for any finite N. Recent numerical studies [32-34] are 
capable, in principle, to bring these large N qualitative 
arguments into a solid theoretical framework. 
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